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Abstract 

This paper addresses new results on the factorization of the general Heun's operator, 
extending the investigations performed in previous works [Applied Mathematics and 
Computation 141 (2003), 177 - 184 and 189 (2007), 816 - 820]. Both generalized 
Heun and generalized Lame equations are considered. 
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1 Introduction 



Fuchsian differential's equations of second order and their confluent forms, 
built for instance from the Ince techniques |I| play a major role in many 
partial differential equations of mathematical physics as Laplace, Helmohltz 
or Schrodinger's equations. 

Without loss of generality the k + 3 regular single singularities can be located 
atx = 0,x = l,x = aj, i = 1, k and x = oo and after an appropriate change 
of function one of the indices at each finite singulary can be shifted to zero 
allowing to write the Fuchsian equation as [2|3] : 
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y (x) + - + + 2^ h/ (a;) + — . = 0. 

\x x - 1 ~{x-ai) x{x - 1) rii=i(^ - di) 

(1) 

The two indices at each singularity are (0, 1 — e$) and (0, 1 — 7), (0, 1 — 5), 
(a, (3) at x — 0, x — 1, and a; = 00, respectively, taking into account the 
Fuchsian relation: 

k 

<* + /? + l = 7 + 5 + J>. (2) 

i=l 



€i = 0, (z = 1, k), gives the hypergeometric equation and k = 1 the Heun's 
equation [7J. The character reducible or irreducible of linear ordinary differ- 
ential equation (O.D.E) is important in relation with the monodromy group 
[2]. We call reducible a second order O.D.E which admits a non trivial solu- 
tion satisfying a first order linear equation, and otherwise irreducible. This 
definition is used by many authors [HI5H6] with or without supplementary con- 
ditions. It allows to give an easy criterion of reducibility of an equation with 
polynomial coefficients, from a factorized approach in the following way. 

Rewriting ([!]) in a polynomial form with D = -p, 



H k [y(x)]= Q k+2 (x)D 2 + Q k+1 D + Q h 



y(x) 



[L{x)D + M{x)\ L{x)D + M{x) 



y{x) 



(3) 



where Qj are polynomials of degree j, possible (or not possible) identification 
of the polynomials L, L, M, M with polynomials Qk+2, Qk+i, Qk will give all 
possible cases of reducibility and, a contrario, criteria of irreducibility for the 
equation TL k [y(x)} = 0. 

This was recently done for the Heun's equation [8|9] 



y"(x) + 



7 



+ 



+ 



X X 



x — a 



y'(x) + 



a/3x — q 



x(x — l)(x — a) 



y{x) 



(4) 



and also for the four confluent Heun's equations [10] ■ For the Heun's equation, 
this factorization is possible and generate 6 non trivial situations in accordance 
with the known F— homotopic transformations [11]. R. S Maier [12] noticed 
that for the Lame equation (k = 1 with 7 = 5 = e = 1/2) 

.. / . 1/1 1 1 \ x l(l + l)x + 4q , , 

2 \i x — 1 x — aj 4x(x — l)(x — a) 



2 



the factorization reduces only to the two cases I = 1, I = ~. He also mentioned 
that for the general Lame equation with k + 3 singularities, the factorization 
is not known. 



Even if some algorithms exist in order to factorize any ODE with rational func- 
tion coefficients ^3ll4|15|fl6] . factorization of an arbitrary k— Heun equation 
is probably untractable, except for using numerical tools. 



The aim of this work is to investigate the case k = 2 for both generalized 
Heun and generalized Lame equations [17J which are also important in lattice 
statistics |18j . These equations have invested a large number of applications 
in physics. See [19J (and references therein) for a nice review on physical 
applications of these equations. To mention a few, retrieved from the indicated 
works, where for instance Heun's equations as well as their solutions impose 
their usefulness, it is worthy of attention to outline their importance in the 
description of quasi-modes of near extremal black branes, in the hyperspherical 
harmonics with applications in three-body systems, in the elaboration of a 
method of calculation of propagators for the case of a massive spin 3/2 field 
for arbitrary space-time dimensions and mass, in parametric resonance after 
inflation, as well as in the separation of variables for the Schrodinger equation 
in a large number of problems, typically for the radial coordinate, and in non- 
linear formulations involving Painleve type equations. A number of traditional 
equations of mathematical physics, as for instance the Lame, spheroidal wave, 
and Mathieu equations, are also particular cases of Heun equations. 



Expansion of the equation (131) gives: 



H k [y(x)} = (LD + M) (lD + M) y(x) 

= [LID 2 + \L (1' + M) + Ml] D + (LM' + MM) } y[ 



x) 



(6) 



and generate the 3 basic relations: 



Qk+2 
Qk+1 
Qk 



LL, 

L (1' + M) + ML, 
LM' + MM. 



(7) 



Computation and properties of the polynomials M = M(x), M = M(x), L = 
L(x), L = L(x) are explicitly given in several situations. 
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2 Preliminary remarks 



1) Qk+2{%) contains k + 2 linear and different factors which are all present in 
the product L{x)L(x) in several ways. For instance, if j factors appear in 
L and k + 2 — j in L, the number of decomposition is given by the classical 
binomial result: 

k+2 

E (D = 2 " +2 ' (8) 
i=o 

excluding the extremal cases 

L{x) = Q k+2 (x), 1 = 1; Q k = 0, 
H k = LD(D + M) (9) 

and 

= 1, L = Qfc + 2(a;), M = (integrable caseM = Qfc+i — Q'k+21 
Q k = M'), H k = D(LD + M) (10) 

which are irrelevant or trivial in this study. 

Selecting a pair L, L inside the 2 k+2 — 2 possibilities, computation of the 
polynomials M(x) and M(x) of degrees m and m, respectively, is done in 2 
steps. 

i) The coefficients in Qk+i{x) and Qk+2{%) allow to linearly compute the 
coefficients of M(x) and M(x), of number m + rh < k from the second 
and third relations in (J7j). 

ii) The third relation in ([7]) gives the value of the accessory parameters 
Pi, {i = 1, k) and the product a/3, using the polynomials M and M com- 
puted in the first step. 

2) The general k— Lame's equation is a k— Heun's equation with 7 = 5 = e = 
1/2, (i = l,k). The indices at each finite singularity are equal to (0, 1/2) 
and the Fuchsian relation ([2]) gives a + (3 = k/2. This equation takes the 
form 



£k[y] 

with now 



Qwy" + Qk+iy' + \ (af3x k + J2 P^^j y = (n) 



1 1 r _-,/ 

Qk+i = -(Qfe+2) =r j (12) 



where 



2" T/ 2 



k 



Q k+ 2 = x(x - 1) - cij) = LL. (13) 
j'=i 
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The property (fl2l) and second relation in (J7|) easily give: 



1 (LV _ LM + ML 

2 VI/ ~ I2 • 



(14) 



The number of situations to be investigated in the decomposition of LL 
can be simplified from some symmetry between polynomials M, M and 
M*, M* corresponding to the decomposition LL, permuting L and L giving 
the factorization: 



AM 



LD + M* 



LD + M* 



y(x) 



A similar relation to (PHI) is 



1 (L\ _ LM* + M*L 

2 \L ~ T 2 ' 



(15) 



(16) 



We deduce from (|T4"|) and ffTBI) that 

L(M + M*) + L(M + M*) = 0. (17) 
Equality ( JTTjl is satisfied if, for example, we have M = — M*, M = —M*. 



3 Lame's equation, (A; = 1) 



Let us consider equation (TJJ with 7 = 5 = e = |. Even if this case has been 
already solved in [7], as a subcase of the full Heun factorization, we easily 
recover the polynomials M and M from the procedure described in section [2] 
for the factorization L = x{x— 1), L example. Unknown polynomials 

M and M must be as: 



M = Ax + B, M = A 



with 



Q 3 (rr) =LL = x 3 - x 2 (a + 1) + xa, (19) 

Q 2 {x) = L(L' + M) + ML = -x 2 - x(a + 1) + | (20) 

Q 1 (x)=MM. (21) 

We easily conclude that 



5 




2 



2 



1 



(22) 



Next the third relation in (J7|) gives: 




(23) 



and the Fuschian relation yields 




2' 



(24) 



From section 2, M* = — M and M* = — M already solve the factorization 
problem with L = x — a, L = x(x — 1). 

A surprising result, consequence of a remark of Maier |12j . is that for L = 
x, or (x — 1), or (x — a), fZ = (rr — l)(x — a), x(x — a), x(x — 1)), the polyno- 
mials M are the same. From the symmetry property M* = —M, M* = —M, 
this result also works for the 3 other factorizations: L — (x — 1) (x — a) , x(x — 
a), x(x — 1); L — x, (x — 1), (x — a). All these results can be seen in table 1, 
presented in the same way as in Ref [9]. 



4 Factorization of Lame's equation /^[y], k = 2 

Equation (pQ) with k = 2, a\ = a, a 2 = b and 7 = 5 = e\ = e 2 = \ is the Lame 



Qi{x) = x(x — l)(x — a)(x — b) 

= x 4 — a; 3 (l + a + b) + x 2 (a + b + ab) — abx, 
< Q 3 ( x ) = \Q' A {x) = \ [4x 3 - 3x 2 (l + a + b) + 2x(a + b + aft) - aft] , (27) 
(^(aO = a/3x 2 + pix + p 2 , 
a + (3 = 1. 

With a first choice L = x{x — 1), L = {x — a) (x — b) , M and M are first degree 
polynomials 



equation 



C 2 [y] = Qi{x)y" + Q 3 (x)y' + Q 2 (x)y = 0, 



(26) 



where 
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Table 1. Factorization of the simple Heun's operator Hi = £i[y] 



x(x — l)(x — a)D 2 + [7(2; — l)(x — a) + 5x(x — a) + ex(x — 1)] D + (a(3x — q)Id 
(LD + M)(LD + M), 7 = S = e = ^. 



L 


L 


M 


M 


a 


/? 


5 


X 


(x — l)(x — a) 


1 

2 


-z + 2±i 


3 
2 


-1 




x-l 


x(x — a) 


1 

2 


-x+§ 


3 
2 


-1 


-(1 + f) 


x — a 


x(x — 1) 


1 

2 


-x+i 


3 
2 


-1 




(x — l)(x — a) 


X 


„, 0+1 
2 


1 

2 


1 
2 


1 




x(x — a) 


x-l 


„ a 
2 


1 

2 


1 
2 


1 




x(x — 1) 


x — a 


X — g 


1 

2 


1 

2 


1 


1 

4 



M = Ax + B, M = Ax + B. 



(28) 



The four equations satisfied by A, B, A, B are given by identification of second 
relation in ([7j) with the second relation in (j2~T|) . The solutions are 



a = i, a = -i, ■* = -!,* = — ■ 

2' ' 2 



(29) 



Comparison of third relation in (J7J) with the third relation in ( 1271) yields 

a(3 = -2 (30) 

and also p\ and p 2 ; 



a + b + 3 
Pi = ^ > 



p 2 



2 

a + & 



(31) 
(32) 



These results as well as those corresponding to the other choices of L and L are 
reproduced in tables 2 to 4 (14 situations) for the cases 7 = 5 = ei = e 2 = ^ 
as in the case k — 1, in which appear essentially 2 situations [12]; the case 
k = 2 generates only 3 situations with a/3 = —2, —3/4, —15/4. 



5 Factorization of the generalized Heun's equation k = 2 



The first extension of Heun's equation with k = 2 can be written as = y) 



y" + 



7 5 
- + 



ei £2 



x x — 1 a; — a x — 6 



x(x — l)(x — a)(x — b) 



with the parameters linked by the Fuchs's condition a + (3+1 = 7 + 5 + ei + e 2 
or in polynomial way 



H 2 [y] = Q 4 (x)y" + Q 3 {x)y' + Q 2 (x)y = 0, 



(34) 



where 



Qi(x) —x(x — l)(x — a)(x — b) 

= x 4 - x 3 (l + a + b) + x 2 (a + b + ab) - xab; 



(35) 
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Qz(x) — j(x — l)(x — a)(x — b) + 5x(x — a)(x — b); 

+ e\x{x — l)(x — b) + e 2 x(x — l)(x — a) (36) 
Q 2 (x) =a(3x 2 + pxx + p 2 . (37) 

The factorized form is 



H 2 [y) = (LD + M) [LD + MJ 

= LLD 2 + (Li! + LM + ML) D + (lM' + MM] . (3£ 



Let us consider a peculiar simple situation 

L = x(x — 1), L = (x — a)(x — b). (39) 

The third relation in $7§: LM'+MM = Q 2 (x) obliges M and M to be of degree 
1. Let M = xA + B, M = xA + B. The unknown A, B, A, B are given using 
also the second relation in (JTj). From Q 2 (x) = Ax(x — 1) + (xA + B)(xA + B), 
we get 

A + AA = ad, 

_ (40) 
-A + _Bv4 + £M = px. 



As in the Lame's case, some kind of symmetry appears when permuting L 
and L in the factorization, allowing to reduce by a factor 2 the number of 
factorizations. Use of the adjoint operator H* of H is now useful. Indeed, the 
Lagrange adjoint Ti,* of Ji is given by [Tf2"Tj 



H* = Q* k+2 D 2 + Ql +1 D + Qt (41) 

with 

= Qk+2, Ql+i = 2Q'k+2 ~ Qk+i, Qt — Q"k+ 2 — Q'k+i + (42) 

It is not very easy to write in general Tt* in the form of the equation (TTTT) 
with star parameters a*, (3*, p*, because the corresponding equations are 
quadratic. But from the factorized form of Ti, computation of polynomials 
M* and M* are trivial from the relations: 

H = (LD + M)(LD + M), H* = (LD + M)*(LD + Mf (43) 
and, from equation fT4"2"|) . we get 
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H* = {-LD + U - M)(-LD + L'-M) 

= (LD + M*) (LD + M*) (44) 



with 



M* = M — L', U' M - /.' (45) 
generalizing the relation (IT51) . 

In tables 5 to 7, we give, in the 14 situations, the values of M, M, a, /3, p\ 
and P2 as well as those corresponding the other choices of L and L. Of course, 
the symmetry in the factorization of TH.2 is not so rich as in the case of £2- 
Nevertheless, a quick sight in the tables allows to classify families having the 
same representation for p\ and p2 in terms of A, A, B, B. See tables 5-7. 
The symmetry between L and L can also be observed in all tables. 



6 Concluding remarks 



Some important features deserve to be pointed out. 

(1) In each case both solutions yi(x), (i = 1,2) of the factorized equations 
are easily obtained solving first the equation 

Ly[(x) + M Vl (x) = 0. 

In all Lame's cases, these two solutions are given by 

y x {x) = y/L(xj, y 2 (x) = Jl{x) [ f X _ . (49) 

J L(x)\J L(x)L(x) 

We recover, of course, for k — 1, the 8 "pseudo polynomials" of Lame 



y{x) = \jL(x)P N (x), (first type), (50) 

Pn{x) being the corresponding Lame's polynomials but now with appro- 
priate parameters a, f3 and pi. 
(2) This approach does not give a new help in the search of a pecular solu- 
tion of a generalized Heun's equation. But it gives a global approach of 
all solutions for arbitrary parameters (e^, Tj). This factorization defines 
in some way a parameter space splitted into many domains defined by 
"factor parameters" given in the several tables. 
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(3) The factorization becomes more and more difficult to handle with the 
increase of the number of k + 2 singularities. Even with k — 3, there are 
already 32 different situations to investigate, what is irrelevant to be pre- 
sented as it implies very long and complicated analytical expressions for 
the polynomials L, M, L, M. A Maple code is now under investigation 
for symbolic computation for arbitrary number k of singularities. 
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Table 2. Factorization of the Lame's operator £2(2/] 
£ 2 [y] = Q4(x)y" + Q 3 (aOl/ + (a/3x 2 + Pl x + p 2 )y = = (LD + M)(LL> + M)y. M = x + B, M = -x + B. 



L 


L 


B 


B 


a 





Pi 


P2 


x(x — 1) 


(x — a)(x — b) 


1 

2 


\{a + b) 


-1 


2 


i(a + 6 + 3) 


~\{a + b) 


(x — a)(x — b) 


x(x — 1) 


-\{a + b) 


1 

2 


2 


-1 


±(3a + 36 + 1) 


-\{a + b + ab) 


x(x — a) 


(x-l)(x-b) 


-\a 


\{b + l) 


-1 


2 


±(3a + 6+l) 


-\{a + ab) 


(x-l)(x-b) 


x(x — a) 


-\{b + l) 


\a 


2 


-1 


i(a + 36 + 3) 


-\{a + Ab + ab) 


x(x — b) 


(x — l)(x — a) 


-\b 


i(a + l) 


-1 


2 


i(a + 36+l) 


-\{b + ab) 


(x — l)(x — a) 


x(x — b) 


-i(a + l) 


\b 


2 


-1 


i(3a + 6 + 3) 


-±(4a + 6 + a&) 



Table 3. Factorization of the Lame's operator £>2[y] (continuation). 
C 2 [y]=Q 4 (x)y" + Q 3 {x)y' + (af3x 2 + Pl x + p 2 )y = = (LD + M)(LD + M)y. M = h 2 + Bx + C, M = ~, a = l, f3 = ~. 



L 


L 


B 


C 


Pi 


92 


(x — l)(x — a)(x — 6) 


X 


-(a + 6 + 1) 


\{a + b + ab) 


±(a + 6 + l) 


-\{a + b + ab) 


x(x — a)(x — b) 


x — 1 


-(a + b) 


\ab 


\{a + b) 


-\ab 


x(x — l)(x — b) 


x — a 


-(6+1) 


\b 




~\b 


x(x — a)(x — 1) 


x — b 


-(o + l) 


\a 


\(a + l) 





Table 4. Factorization of the Lame's operator C 2 [y\ (end). 
C 2 [y] = Q 4 (x)y" + Q 3 (x)y' + (af3x 2 + Pl x + p 2 )y = = (LD + M){LD + M)y. M = ~x 2 + Bx + C, M = 1 a = ~, = | 



L 


L 


B 


C 


Pi 


P2 


X 


(x — l)(x — a)(x — b) 


a + b+l 


-\{a + b + ab) 


|(a + 6 + l) 


-i(a + 6 + a6) 


x - 1 


x(x — a)(x — b) 


a + b 


-\ab 


|(2 + a + 6) 


-i(4a + 46 + a6) 


x — a 


x(x — l)(x — b) 


b+l 


-\b 


§(2a + 6 + l) 


-i(4a + 6 + 4a6) 


x — b 


x(x — a)(x — 1) 


a + 1 


-\a 


|(a + 26 + l) 


-±(a + 46 + 4a&) 



Table 5. Factorization of the Heun's operator 7Y 2 [y] 
n 2 [y]=Q±{x)y" + Q 3 (x)y' + {af3x 2 + Pl x + p 2 )y = = (LD + M)(LD + M)y. M = Ax + B, M = Ax + B. 



L 


L 


A 


A 




B 


a 


/3 


x(x — 1) 


(x — a)(x — b) 


5 + 7 


-2 + ei + e 2 


-7 


a — e 2 a + b — e\b 


-2 + ei + e 2 


7 + <5 + l 


(x — a)(x — b) 


x(x — 1) 


ei + e 2 


-2 + 7 + 5 


— e 2 a — ei6 


1-7 


ei + e 2 + 1 


-2 + 7 + 


x(x — a) 


(x-l)(x-b) 


ei + 7 


5 - 2 + e 2 


—7a 


1 - e 2 - 5b + 6 


5 - 2 + e 2 


7 + ei + 1 


(x - l)(x - b) 


x{x — a) 


5 + e 2 


-2 + ei +7 


-56 - e 2 


a — 7a 


1 + <5 + e 2 


-2 + ei + 7 


x(x — b) 


(x — l)(x — a) 


7 + ^2 


5-2 + ei 


—76 


1 — ei — 5a + a 


£ -2 + ei 


1 + 7 + e 2 


(x — l)(x — a) 


x(x — b) 


ei + 5 


e 2 - 2 + 7 


—5a — e\ 


6 — 76 


1 + ei + (5 


e 2 - 2 + 7 





L 


L 


Pi 


P2 






x(x — 1) 


(x — a)(x — b) 


-A + BA + AB 


BB 




(x — a)(x — 6) 


x(x — 1) 


-{a + b)A + BA + AB 


abA + BB 


x(x — a) 


(x - l)(x - b) 


-aA + BA + AB 


BB 


(x - l)(x - b) 


x(x — a) 


-{l + b)A + BA + AB 


bA + BB 


x(x — 6) 


(x — l)(x — a) 


-bA + BA + AB 


BB 


(x — l)(x — a) 


x(x — 6) 


-(1 + a)A + BA + AB 


aA + BB 



Table 6. Factorization of the Heun's operator 7Y 2 [y] (continuation). 
n 2 [y] = Qi{x)y" + Q 3 (x)y' + (a(3x 2 + Pl x + p 2 )y = = (LD + M){lD + M)y. M = Ax 2 + Bx + C, M = A. (47) 



L 


L 


A 


5 


C 




a 





Pi 


P2 


(x — l)(x — a)(x — b) 


X 


ei + 5 + e 2 


—5a — 5b — e 2 — e 2 a — ei — ei& 


e 2 a + ei& + oafe 


-1 + 7 


ei + 5 + e 2 


-1 + 7 


BA 


CA 


x{x — a){x — b) 


x - 1 


ei + 7 + £2 


— e 2 a — 7a — ei& — 76 


7a6 


-l + o 


ei + 7 + e 2 


-1 + 5 


BA 


CA 


x{x — l)(x — b) 


x — a 


<5 + 7 + e 2 


_ e2 _ ry _ ^5 _ J5 


76 


-1 + ei 


<5 + 7 + e 2 


-1 + ei 


bA 


CA 


x{x — a){x — 1) 


x — b 


ei + 5 + 7 


—7a — 5a — 7 — ei 


7a 


-l + e 2 


ei + 5 + 7 


-l + e 2 


BA 


CA 



Table 7. Factorization of the Heun's operator 7Y 2 [y] (end). 



n 2 [y] = Qi{x)y" + Q 3 {x)y' + (a(3x 2 + Pl x + p 2 )y = = (LD + M)(LL> + M)y. M = Ax 2 + Bx + C, M = A. (48) 



L 


L 


A 


A 




C 


X 


(x — l)(x — a)(x — b) 


7 


-3 + ei + (5 + e 2 


2a + 2b + 2 — 5a — 5b — e 2 — e 2 a — ei — e\b 


—ab — a — 6 + e 2 a + ei6 + Sab 


x-1 


x{x — a)(x — b) 


5 


7 - 3 + ei + e 2 


2b — e\b — 7& + 2a — 7a — e 2 a 


—ab + 706 


x — a 


x{x — l)(x — b) 


ei 


-3 + <5 + 7 + e 2 


26 - (56 - 76 + 2 - 7 - e 2 


-6 + 76 


x — b 


x{x — a)(x — 1) 


e2 


-3 + 61+7 + 5 


2 — ei — 7 + 2a — 7a — 5a 


—a + 7a 



L 


L 


a 


/? 


Pi 


P2 


X 


(x — l)(x — a)(x — 6) 


-3 + ei + (5 + e 2 


2 + 7 


B + AB 


AC 


x — 1 


x(x — a)(x — 6) 


7 - 3 + ei + e 2 


2 + 5 


-2A + B + AB 


-B + AC 


x — a 


x(x — l)(x — 6) 


-3 + 5 + 7 + e 2 


2 + ei 


-2aA + B + AB 


-aB + AC 


x — 6 


x(x — a)(x — 1) 


-3 + ei + 7 + (5 


2 + e 2 


-2bA + B + AB 


-bB + AC 



